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Abstract. Using field-theoretic methods, we calculate the internal energy for the One-Component Plasma 
(OCP). We go beyond the recent calculation by Brilliantov [N. Brilliantov, Contrib. Plasma Phys. 38, 
489 (1998)] by including non-Gaussian terms. We show that, for the whole range of the plasma parameter 
r, the effect of the higher-order terms is small and that the final result is not improved relative to the 
Gaussian theory when compared to simulations. 

PACS. 52.25.Kn Thermodynamics of plasmas - 61.20.Gy Theory and models of liquid structure - 05.20.-y 
Statistical mechanics 

1 Introduction matter physics, jellium is often used as a reference state 

when calculating the electronic structure of solids. When 



In its simplest form, the One-Component Plasma (OCP) 
I 

'■^ consists of a collection of n equally charged point-particles 

o 
o 



immersed in a neutralizing background that assures global 
charge neutrality of the system. The OCP is important 



generalized to a Two-Component Plasma (or Restricted 
Primitive Model, if hard-core interactions are taken into 
account) , it can describe electrolytes and electrostatically 
stabilized colloidal solutions. For reviews see [|l|,D. 



[ in several areas of physics as a starting point from where 

■ " " ' concepts or more refined theories are derived. For instance. Different analytical techniques were employed in order 

in astrophysics both OCP and its quantum-mechanical to understand the OCP. These were, in most cases, based 

counterpart (the electron gas or "jeUium") are used in on integral equations (as for instance in ||,^,^) or mod- 

the description of degenerate stellar matter (interior of ified Mayer expansions i.e., low density expansions 

white dwarfs and outer layer of neutron stars) and the that used infinite rcsummation of diagrams that accounted 

interior of massive planets Hke Jupiter^]. In condensed for the long range character of the Coulomb interaction, 

e-mail: amoreira@mpikg-teltow.mpg.de Comparison of the theoretical results with experimental 
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data is usually not possible. It is here where simulations where A is an arbitrary constant, v{r) = is/r is the bare 

(or "computer experiments" 1^) play a particularly impor- Coulomb operator and £b = e^/4 tt e Kb T is the Bjerrum 

tant role, by providing a test ground where the suitabil- length (the length at which two elementary charges have 

ity and range of validity of the different approaches can an interaction energy equal to the thermal energy). The 

be checked. In general, the simulations use Monte Carlo charge density Pc{y) is defined as 

technique P, |9|, p^, pi] and yield quantities like the internal " 

p,(r) = q^(5(r-r,), (2) 

energy or the pair distribution function g{r). j=i 

In a recent paper ||l^, a field-theoretic approach was where q is the valency of the particles and p_ is the 

used to treat the OCP. The introduction of a cut-off at uniform density of the neutralizing background. Global 

small wavelength (large-fc) related to the mean distance charge neutrality implies that /9_ = qn/V , where V is the 

between particles led to a good agreement, for all values of volume of the system. The Coulomb self-interaction of the 

the plasma parameter between the calculated internal particles is given by 

energy and simulation results. The field-theoretic action 2 n"^ r, r A tt P 

IP q n q n I qk 'iirig 

Escif = ^ viO) = — / -—3 — ^. (3) 

used in neglects terms other than the Gaussian ones, J \ ) 

i.e., it goes up to second order in the fluctuating field. This corresponds to an infinite shift in the chemical poten- 
Here we extend this by including more terms in the action tial, which is unimportant for the thermodynamics of the 
and calculating consistently, using the same cut-off, their system. We keep this term for reasons that will become 



contribution to the internal energy. As our main result, we clear later (cf. Eq. 10 below). 

show that the higher order terms do not affect significantly We can apply the Hubbard-Stratonovich transforma- 
the results obtained with the Gaussian theory derived in tion and obtain a partition function that depends on a 
. fluctuating fleld <j)- This follows closely what has been done 

in and so we state here the final expression 

2 The Field-Theoretic Model „ r 06 ( 

Let us assume a classical system where n positively charged — 1 y dr dr' 0(r) u^jlj (r — r') (j){r') + W[(I)]^, (4) 

particles are immersed in a neutralizing negatively charged 

where S = — nln(c A^) is the ideal entropy of the particles 

background. The partition function of this system is 

(c — n/V) and 



ni 



n/ A3 _ ^^Pj^^^if J ^^exp|-i J drdr'(/.(r)z;-i(r-r')</)(r')| 

-iy'drdr'p,(r)«(r-r')pc(r')| (1) 
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We defined Z°'^ as the "excess partition function" , the part as propagator. The term 

of Z that accounts only for the interactions between the -^dh _ ^ ■^dh 

Kb! \ Zq 

particles. The propagator v^^ in (|^) is the inverse of the 



bare Coulomb operator and u^jlj in (Q) is given by 2 tt 

v^Ur)^-^^ + q'c5{v). (6) 



V_ 

2 





1 , 



(10) 



-1 f \ _ V^(5(r) 2 xf \ ^R\ is the Debye-Hiickel contribution to the free energy. Notice 



It is easy to show that this propagator is the inverse of the 
Debye-Hiickel operator WDH(r) — ^se'^'/r, where is 
the screening length given by — AttIb c. The W[(j)] 
is an infinite series in (f> that contains only non-Gaussian 
terms. Up to eighth order it reads 



that £'soif automatically regularizes this integral in the 
ultra-violet, allowing its evaluation without the need of a 
small wave-length (large-fc) cut-off. 

Using (|l|) and neglecting the term log(e^["*l) in (|), 
we get the excess free energy per particle 



3! 4! 

,1/3. 



/g y ^_ where = (7^ i?B (4 TT c/3) is the (diniensionless) plasma 



10 ( 



5! 



parameter. From the excess free energy, we can get the 



6! \ / internal energy per particle 

\IrV 



'-21(/)2,/,5-7O03 04 + 21O^2 ,^3 jj Qj,^ V3 V2 

7' V J u=—-- = r-L- = -\r^'\ (12) 

, hKb 1 oi 2 
(^-28 06 02_5g^^_35042 

V This is the resulting u for what we call from now on the 

+ 420 04 ,/)2% 560 02 - 63002^^ "Gaussian theory without cut-off" . 

(7) Eq. (^ij) is the well known Debye-Hiickel limiting law, 

with = c and i2 = -1; to simplify the notation, we which is asymptotically exact for vanishing T. At T 

use 0" = / dr 0(1) the expression ( p^ ) already yields poor results when 

If we define compared to simulations (cf. Fig. lb). At large F (cf. 



^DH = I Dcj) exp|-i J drdrXr) i;DH(r - r') 0(r')|, 

(8) 

then the excess free energy of the OCP is given by 



1 ^ ^ ^ ^ ^ ^ ^ ^ /N ,/ /\1 ^'^S- Isi) this inadequacy is particularly clear: fits to sim- 
ulation data show a linear behavior ||l|,|l^ in the internal 
energy, and not a 3/2 power law. 

Brilliantov||l2j calculated the Gaussian theory as de- 

^™ — 1 _ p ^ f ^DH \ , / w[4>]\ 

rp ~ ~ ^*^S ^ ^ — -t^scif — iog I I — iog\e j . picted above but introduced a modification, namely a large- 

(9) k cut-off. This is justified with ideas that follow the Debye 

The angular brackets correspond to a Gaussian average theory for the specific heat in solids, stating a direct rela- 

where the the inverse Debye-Hiickel operator (ph is used tion between the number of allowed k modes in the system 
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and the number of degrees of freedom 3 n. The allowed 
wave vectors would be approximately inside a spheres of 

1 /3 

radius ko = (9c7r^) , which is used to substitute the oo 
in the integral in The agreement between the inter- 
nal energy obtained with this cut-off and the results from 
Monte Carlo simulations ]lO| , [Tl| are good (cf. Figs. 1 and 
2). What we will show next is that the inclusion of con- 
tributions up to eighth order in coming from the term 
log(e^t*l) does not change significantly this result. 

When the cut-off ko is used in ([l0| ) , we get the Debye- 
Hiickel excess free energy per particle 



/dh — 



-I (or; 



3/2 



arctan 



hr -\og{i + br)] (13) 



1 /3 

where F is the plasma parameter and h = (Z/tt^) 2/3 
(keeping the notation used in jlj] ) . Notice that in the limit 
r ^ Q this expression reduces to (pi]), as it should. This is 
the excess free energy used in [|2) to calculate the internal 
energy. 

In order to go beyond the Debye-Hiickel level, we do 
the cumulant expansion 



log(e^ 



1 



(14) 



where 



b-') = ^ [dk-^ 







96r 



(16) 



1 — vbr arctan 



/bT 



and 

X„^ = 4 7r / dr (0(0) </)(r))^ 



4 ^9^^ ^™-3(^^)3(™-l)/2 



'r(3 



/6r). 

(17) 



r{m, x) is the incomplete gamma function [p^ and a — 

1 /3 

7r/(9c7r^) is a small distance cut-off. Ideally, the in- 
tegrals Xm should be performed in fc-space with the mo- 
mentum cut-off given by ko. However, for m > 3, the 
Fourier transformed integrals cannot be solved in a closed 
form; on the other hand, the integrals, when written in 
real space are not difficult to calculate exactly, provided a 
small distance cut-off. Since a large-fc cut-off corresponds 
in real space roughly to a small distance cut-off, we intro- 
duced as an approximation the small-r limit a such that 
ko — I a 

Putting (|lj) and ( [l5| ) into (^ , we finally obtain the ex- 
pression for the excess free energy of the OCP with con- 
tributions up to eighth order in 0. The internal energy 
follows as in (|l^). 



Using (0) and going up to eighth order in (/), we obtain 3 Results and DisCUSSion 



log(e 



108 



X3 + t: W / Xi 



^(0^X1-^X4 (15) 



In Figs, la and lb we show u obtained from simulations [pO|, 
|Tl| (black circles), from the Gaussian theory wit/i||l^ and 
without (Eq. p2|) cut-off (respectively dashed line and dash- 
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dotted line) and from the results obtained here with the remaining ones to the eighth order. We can then write 
higher order terms (full line). down u as a sum of three terms, viz. 

In the strongly coupled regime (F > 1, Fig. la), the , a , a /1n^ 

inclusion of higher order terms does not affect the results 

where wdh is the Gaussian contribution to the internal 
energy coming from (p^. In Fig. 3 we plot Z\m6/|udh| 
and zAws/IudhI- As we can see, Auq is approximately one 
order of magnitude smaller than udh and Aug two order 
of magnitude smaller than udh • We expect that the inclu- 



obtained with the Gaussian theory with cut-off; both re- 
sults are indistinguishable on this scale. We also calculate 
the relative error in u, defined as 

u — Us 

RciT — — i — [— , (18) 



where the subscript s stands for simulation. In Fig. 2a we sion of terms of order higher than eight will not change 

show, for r > I, i?crr for the Gaussian theory with cut-off significantly the picture given here. 

(black circles) and for u with the higher order terms calcu- In summary, we have calculated higher order contribu- 

lated here (white circles). The agreement between theory tions to the internal energy of the OCP. We have shown 

and simulation is good, with deviations between —2% and that (i) the effects of these higher order terms are small 

2%[|l^. Notice however that for F < 10 the inclusion of relative to the previously calculated Gaussian theory [p^ 

the higher order terms make the results worse, relative and (ii) they do not improve the agreement between the- 

to the Gaussian theory with cut-off, when compared to ory and simulation This shows that the Gaussian theory 

simulation. with the cut-off kg introduced such as to approximately 

In the weakly coupled regime (-T < 1, Fig. lb), this include strong nearest-neighbor correlations in the high- _r 

trend is confirmed. In Fig. 2b it is clear that the devia- limit is very accurate for describing the 0CP[p7|. 
tion between theory and simulations are larger when the Our calculation consists of two major steps, viz., the 

higher orders are included. However, in both cases i?orr in- expansion of the excess free energy (^ in cumulants of 

creases with decreasing F. This is not surprising though, and the introduction of the cut-off ko suggested in Ref. 

since u goes to zero as F decreases, making the relative er- [ p^ . In principle, the first step can be improved system- 

ror very sensitive to small differences between theory and atically by including higher order terms, while there is no 

simulation. clear recipe for improving the second one. The results we 

In Fig. 3 we assess the importance of the higher order obtained here are then a consequence of the approximate 
corrections computed here in comparison to the Gaussian way of calculating the cut-off. Since the precise value of 
theory with cut-off. The higher order terms of the excess kg is not uniquely determined, it may be treated as a fit 
free energy are given by (p^): the first term in the rhs parameter. By fitting the theory to simulation at the high- 
corresponds to the sixth order in (j) correction and two F limit we obtain k^t — 4.417 c^/'^, which is close to the 
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value ko — (9 tt^ c) ~ 4.462 c^/^ used by Brilliantov and 16. Although the deviations are small, we should notice that 



also used by us. 
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they are statistically significant since the numerical er- 
ror given by simulations in this range of F is below 
±0.02%g. 

17. As a side remark, we note that the application of the cut- 
off ko to the Two-Component Plasma (TCP) leads to un- 
physical results. At the Debye-Hiickel level, the excess free 
energy per volume goes for small c = (n+ -I- n-)/V as 
f™ —[4-K£Bq'^c)^^^/12n (the already discussed Debey- 

Hiickel limiting law) and as f''^ (9/7r)^/3 ^ijg^c^/a for 

large c; at this limit total free energy per volume is un- 
bounded from below and non-convex. At finite tempera- 
tures, the Maxwell construction leads to the coexistence 
between states with a zero and an infinite density, which 
is clearly unphysical and signals that this approach is not 
appropriate for the TCP. 

This difference between the OCP and the TCP becomes 
more transparent in a recent systematic low-density expan- 
sion of the free energy, where the TCP expansion coeffi- 
cients are shown to be dominated by ion-pair formation 

18. R. R. Netz and H. Orland, to be published. 
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Fig. 1. Internal energy (u) as function of the plasma parameter 
r in the range (a) 1 < T < 300 and (b) < T < 1. The 
full line denotes the result obtained here, the dashed line the 
Gaussian theory with cut-off jl2j, the dash-dotted line the u 
obtained from the Debye-Hiickel limiting law (Eq. |l^). The 
points denote simulation results from Refs. jioj (0 < -T < 1) 
and [0 (1< r < 300). 

Fig. 2. The relative error _Rcrr (Eq. |l^) as function of F for (a) 
1 < r < 300 and (b) 0.1 < T < 1. The white circles represent 
the relative error for the u obtained here with the higher order 
corrections; the black circles represent the relative error for the 
Gaussian theory with cut-off 

Fig. 3. Comparison of the sixth order and eighth order cor- 
rections with itDH- Notice that Aue is roughly one order of 
magnitude smaller than udh and Auq is two orders smaller 
than itDH- 



Moreira and Netz, Fig. la 




Moreira and Netz, Fig. lb 




Moreira and Netz, Fig. 2a 




Moreira and Netz, Fig. 2b 




Moreira and Netz, Fig. 3 




